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MINIMUM-WAVE-DRAG AIRFOIL SECTIONS FOR ARROW WINGS 
By Morton Cooper and Frederick C. Grant 

SUMMARY 


A linearized theoretical, analysis has been made to determine 
mlnlmum-wave-drag airfoil sections for arrow wings having the same air- 
foil sections at all spanwlse stations. The drag of the wings was mini- 
mized subject to the condition of either a given thickness ratio at a 
specified chordwise location or a given vol'ume. Numerical computations 
of the airfoil shape and wing wave drag were performed for a delta wing 
and for an arrow wing having a ratio of the tangent of the tralling-edge 
sweep angle to the tangent of the leading-edge sweep angle of 0.4. The 
range of the ratio of the tangent of the leading-edge sweep angle to the 
tangent of the Mach line sweep angle extended from 0 (two-dimensional) 
to 2 . 5 . 


INTRODUCTION 


At supersonic speeds^ the theoretical determination of minimum- 
wave-drag wings is of importance in establishing criteria for efficient 
aeronautical design. There is a rapidly Increasing amount of information 
based on linear theory (refs. 1 to 4) and slender-body theory (refs. 5 
and 6) to serve as a guide for the selection of three-dimensional 
minlmmn-wave-drag wings of given thickness or volume. As yet, however, 
little information exists for specifying the shapes for minimum-drag 
airfoil sections for wing plan forms of current interest at supersonic 
Mach numbers. It is the p\nrpose of the present paper to study this 
problem by linearized theory. 

A constant-thickness-ratlo wing of arrow plan form with an arbitrary 
airfoil section is considered. For convenience, the root chord of the 
wing is assumed to be \inity. The airfoil section (assumed to be the 
same at a.i l spanwlse stations) is composed of a series of strai^t-llne 
segments. The drag of such a wing is then minimized subject to the 
restriction of either a given thickness ratio at a specified chordwise 
location or a given voliome, and the equations for the corresponding 
minimum-drag airfoil shapes are determined. Computations are made to 
establish the minimum number of sides required for a reasonable deter- 
mination of the airfoil shape and drag. Optimum sections and the 
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corresponding wing drag coefficients ajre then calculated for a series of 
delta and arrow wings swept various amounts ahead of and behind the Mach 
lines. As a limiting case^ the results of the optimum-airfoil calcula- 
tions for a delta wing of fixed volume are compared with those of the 
corresponding airfoil derived analytically from slender-hody theory . 

The authors wish to acknowledge the cooperative efforts of the 
National Bureau of Standards personnel, particularly Otto Steiner, 

Sally Tsingou, and Norman Levine, who set up the coding and carried 
through the numerical computations on the National Bureau of Standards 
Eastern Automatic Computer during June 1955 » 


SYMBOLS 


M 


P 

t 

N 

k 


ki 

%+l 
n = k/p 

x,y,z 

A 


Amax 

A 

z(x) 

K 


Mach number 

cotangent of Mach angle, /m^ - 1 ... .... 

thickness ratio 

number of equally spaced straight-line segments used to form 

one side of symmetrical airfoil (number of chord divisions) v 

tangent of sweep angle of any arbitrary ridge line 

tangent of sweep angle of leading edge 

tangent of sweep angle of trailing edge 

Cartesian coordinates (see fig. l) 
cross-sectional area of body at station x 
maximum cross-sectional area of body 
slope of airfoil surface divided by thickness ratio 
airfoil shape function 

ratio, of airfoil cross-sectional area to area of circum- 
scribing rectangle 


P 


Integer specifying last subsonic ridge line 


V 
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m integer specifying location of maximum thickness 

l,j,r,a ;|3 arbitrary indices 
AAj* = “ ^r~l 

B' functions used to determine drag 

F function of AK defined by equation ( 8 ) 

Cp drag coefficient based on ving area 


ANALYSIS 
Linear Theory 


In reference 1 , the basic equations for the drag of arrow wings of 
double-wedge airfoil sections have been obtained by the superposition 
of constant-strength source distributions within the wing plan form. 
This method has been extended in reference 7 into a generalized proce- 
dure whereby the wave drag of arrow wings having arbitrary profiles may 
be determined. This extension is accomplished by using a finite number 
of constant-strength source distributions and hence entails approxi- 
mating the airfoil section by a multisided polygon. If the sides of 
the polygon are equally spaced (fig. l) and the airfoil section is 
symmetrical., the drag at zero angle of attack is, in sll^tly different 
notation from that of reference 




e(. - 


- ± ± Cj - i)‘ 

i=l i=i 




AAi AAj - 


n 


p N+1 
N N +1 

r F 

i=p+l 


V 


( 1 ) 
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ji 

■where the new variable AA3- is given by 

* 


and where A is defined hy the relationship 


dz 



In equation (l), the index p defines the last subsonic ridge line (a 
subsonic line is one whose normal ffech number component is subsonic; 
see fig. 1) and has the integer value determined from 


- 1) , j ' j. > 

^1 " ^Wfl ^1 - ^K+1 


( 2 ) 


The values of n are geometrically related by the equation 


ni = ni - 


N 


( 5 ) 


The foiTOulas for the functions b', D, and f obtained from refer- 

ence 7 


B’ 




cosh'^i 


+ 
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In order to Investigate optimum airfoil sections for arrow wings, 
the drag, as given by eqjuation (l), is minimized subject to the speci- 
fied thickness or volume condition and the airfoil shape parameters AXi 
are determined. 
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The specified conditions are: 


(1) The airstream ahead of and behind the wing is parallel to the 
wing- chord plane . . 


01 = 



AXi = 0 


(5) 


(2) The airfoil most close at the leading and trailing edges 

N+1 

{Zi2 = XI ^ ^1 = ^ 

1=1 


(3) The airfoil must have a given thickness ratio t at a specified 
fraction of the choixi m/N 


03a = 2 ^ (m + 1 - i)AAi - N = 0 


(Ta) 


or the airfoil must have a given cross-sectional area so that the wing 
has a given volume 



l^AXi - K = 0 


(71=) 


In equation K specifies the ratio of^the cross-sectional area 
of the optimum airfoil to the area of the circumscribing rectangle (a 
point discussed subsequently) . 


The shape of the minimum-drag airfoil can be obtained analytically 
by the extreme-value theory of reference 8, Since the thickness dis- 
tribution z(x)/t of an optimum airfoil section is independent of the 
thickness ratio (a fact that can be inferred from eq.. (l)), the minimum- 
drag shapes are determined by first formulating the function 


F 




4 - 



+ X01 + Y02 + Z 



( 8 ) 


*The braces Indicate that the appropriate condition, either 03 q_ or 
03^, is used. 
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where 


8 





t2 


is the function of AA specified hy equation (l) 


and 02^ 03a^ 05b functions of AX specified by- 

equations ( 5 ), ( 6 ), ( 7 a), and (7b), respectively. The variables X, Y, 
and Z are undetermined multipliers. The optimum. AX values are then 
obtained by solving the system of H + it linear simultaneous algebraic 
equations given by; 


and either 


or 


=0 

8AXr 

(9) 

0 

II 

H 

(10) 

0 

II 

OJ 

(11) 

05a ~ 

(12a) 

It 

0 

(12b) 


Equations ( 9 ) to (l2) are sho-wn in matrix notation in exploded 
form in figure 2 ; the relations for the individual entries in the matrix 
are also Included. The interchange of the area condition for the thick- 
ness condition aimply changes one row and one column of the matrix as 
indicated in figure 2 . 

In the nimierical evaluation of the optimum airfoil for a specified 
area, the parameter K depends upon the airfoil shape. The calc\ilated 
values of AX, however, are directly proportional to any assumed value 
of K. Hence, any arbitrary value may be initially selected and the 
solution linearly scaled so that the relationship for the airfoil 
coordinates 


*” 2z-p_ ~i — 


t AXj^^ q 
N 


(13) 


yields a value of Zj./t equal to O .5 a’t the location of maxi mum thickness. 
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Slender-Body Theory 

In reference 6 , it is established, on the basis of slender-body 
theory, that the optimum wing and body of revolution (both for a given 
length and volume) have the same axial cross-sectional area distribution. 
Hence, as a limiting case for n arbitrarily large, the airfoil section 
obtained by this slender-body method can be con^iared with the results 
of the previous linear analysis. 

The shape of the airfoil section of a constant-thickness-ratio delta 
wing having the same axial cross-sectional area distribution as a body 
of revolution is 


z(x) ^ 1 
t 2 


(1 - x) ^^— + 2 A(x) 


dx 


/t ,, N 

dA(x) 

“ ^%axy 

_dx 


+ 2A 


(%iax} 


(ih) 


where Xm»Y is the location of maximum thickness of the airfoil and 
a(x) is the axial cross-sectional area distribution of the body or wing. 
Equation U^-), when applied to the minimum-drag body of revolution for a 
given len^h and volume (ref. 9) 


Mil = 8(x - x 2 ) 5 /= 
Afflax 

results in the airfoil shape 


z(x) 

t 


1 

2 


(x - x^)^^^(2x - 3)(x - 1) 

(^•%ax " ^%ax^) (^^^max ■■ ^)(%iax " 


(15) 


(16) 


V 5 - ]/5 

where Xj^g^ = ^ • 

The shape determined by equation (16) is compared with the linear- 
theory shapes in the following section. 
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NUMERICAL RESULTS 


Effect of number of sides N .- The initial problem consisted in 
establishing the number of sides per surface N required for a reason- 
able appro x i ma tion to the shape and drag of the limiting airfoil having 
an infinite niunber of sides. For this trial calculation a delta wing 
was assumed and the thickness condition applied. A location of maximum 
thickness m/N of 0.20 of the chord vas selected and a value of n^ = I .7 
(subsonic leading edge) was assumed. A subsonic-leading-edge condition 
was selected since it was anticipated that the limiting airfoil would 
have a round nose and hence require the largest number of flat sides 
per s-urface N, for a good approximation. 

The resiilts of the computations for N = 5 , 10, 20, and 14-0 sides 
are shown in figure 3 and presented in table I. Perhaps more striking 
in figure 3 than the Influence of the number of sides is the bump 
located at the position of maximum thickness. From an examination of 
all the ninnerical conrputations, it can be inferred that this bump, 
which tends to become more cusplike as the number of sides Increases, 
occiirs only when the thickness auxiliary condition is applied with a 
subsonic maximum- thickness line. For such conditions, the resulting 
airfoil is in violation of the assimptions of linearized theory 
hence represents a restriction on the use of such a theory. Further- 
more, in the use of such airfoils, flow separation would undoubtedly 
be a governing factor. As is demonstrated subsequently, if the maximum- 
thickness line is supersonic or if the area condition is used, no such 
difficulty arises. From figure 3, it can be seen that an Increase in 
the number of sides from 20 to IfO has little effect on either the drag 
or the airfoil shape; a value of N = 20 was therefore selected for 
all succeeding computations. Since in this test case the limiting 
airfoil appears to be a difficult one to approximate by line segments, 
it is believed that the value of N = 20 is entirely adequate.’ 

Range of cocputations .- The conputations of the wave drag air- 
foil shape for the optimum airfoil sections were performed for delta 
wings and for arrow wings, that is, wings having a ratio of the tangent 
of the trailing-edge sweep angle to the tangent of the leading-edge 
sweep angle of 0 and 0.4, respectively. The range of values of nj_ 
extended from 0.5 to 2 . 5 . For the thickness condition, maximum- 
thickness locations of 0.20 and 0.50 of the chord were Investigated. A 
complete summary of the range of conputations is presented in table II. 

The final computations are tabulated in tables III and ly for the delta 
and arrow plan forms, respectively, and are plotted in figures 4 and 5 . 
Figures 4 and 5 also Include the two-dimensional (n^ = O) optimum resiilts. 

Optimum section for the thickness condition .- Figure 4 presents a 
comparison of the optimum- airfoil-shape and drag calculations with 
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corresponding double- vedge results (ref. 10 ) for two chordwlse locations 
of maximum thickness: 0.20 and O. 5 O. For -the 0.20 location of maximum 

thickness (fig. i<-(a)) all results for n^ > 1 correspond to subsonic 
maximum- thickness linesj whereas for the O. 5 O location (fig. 4(b)) all 
results correspond to supersonic maximum- thickness lines. 

For both locations of maximum thickness (figs. 4(a) and 4(b)), plan 
form has little effect on the optimum shapes. For n^ < 1 the optimum 
airfoil has almost the same shape and drag as the double-wedge airfoil 
(figs. 4(a) and 4(b)). For the 0.50 location of maximum thickness 
(fig. 4(b)), both wings have, for subsonic leading edges (n]_ > l) , opti- 
mum sections which are more full at the front, more closely resemble 
subsonic airfoils, and have about 16 percent less drag than the corre- 
sponding double-wedge section. For the 0.20 location of maximum thick- 
ness (fig. 4(a)), the progressive decrease in magnitude of the bump at 
maximum thickness is apparent as n 2 , decreases. No discussion of the 
drag of wings having subsonic maximum-thickness lines is considered to 
be warranted in view of the limitations on the theory and the separation 
problem in a real flow for this case. 

Optimum section for the area condition .- Figure 5 presents the drag 
and optimum airfoil shapes for delta and arrow wings of a given volvime. 
Again the airfoil shapes are plotted at the appropriate n 2 _ values. As 
was the result for the thickness condition, plan form had little effect 
on the optimum section, and for the supersonic leading edges, the two- 
dimensional optimum (the biconvex section) results for a3J. practical 
purposes . This latter result was established for an elliptic-plan-form 
wing in reference 5. For subsonic leading edges (n]_ > l), the airfoils 
become more full at the front and again resemble subsonic airfoils. The 
maximum- thickness location has moved forward for large values of n]_ 
and for the last case computed occurs at 25 percent of the chord. The 
progressive rounding of the nose section of the delta wing for hl^ly 
subsonic leading edges is shown in figure 6 where the slender-body air- 
foil (n 2 _ — ?•“) given by equation (16) (table V) is compared with airfoils 
derived from linear theory. For this slender-body airfoil, the location 
of maximum thickness is at I 9 percent of the chord. For a value of nq 
of 2.5 the linearized result already indicates a close similarity with 
the limiting result (nq-^“) as obtained by slender-body theory. 

The optimum linear- theory airfoil, subject to the thickness condi- 
tion, is not compared with the slender-body airfoil (derived from an 
optimum body of given length and diameter) because the auxiliary condi- 
tions are not comparable. 

Comparison of airfoil sections for wings of given volixtne .- In order 
to establish a basis of comparison for the drag results in figure 5^ ’the 
drag of wings of the same plan form and volume but having biconvex and 
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NACA 65A-Beries sections was calculated (again for N = 20) and is com- 
pared in figure 7 with the drag of optimum sections. This comparison is^ 
perhaps, somewhat artificial in that preselected values of the location 
of maximum thickness exist for the biconvex (50 percent) and the NACA 
65A-series (lj-0 percent) wings. These results are believed to be of 
interest, however, inasmuch as the biconvex section is the two-dimensional 
(ni = 0 ) optimum for a given area and the 65A-series section is of current 
interest . 

For the supersonic leading edges, the drag reduction of the optimum 
section as compared with that for the biconvex section is small a.nf^ 
amoimts to less than 10 percent. For the subsonic leading edges, the 
drag reduction of the optimum section is 2 ^ percent of that of the 
biconvex wing of delta plan form and 33 percent of that of arrow plan 
form for a value of nj^ of I .7 (6o° leading-edge sweep at M = l.lt, for 
exan 5 )le). Corresponding values of the drag reduction of the optimum 
section as compared with those for the NACA 65A-series sections would be 
about 10 percent and I 5 percent, respectively, at the same value of n-j_. 


CONCLCJSIOl© 


A linearized theoretical analysis has been made to determine the 
miniraum-wave-drag airfoil sections for arrow wings having the same 
airfoil sections at all spanwlse stations. The drag of the wings was 
minimized subject to the condition of either a given thickness ratio 
at a specified choidwise location or a given volume. Numerical computa- 
tions were performed for a delta wing and an arrow wing having a ratio 
of the tangent of tralllng-edge sweep angle to the tangent of the leading- 
edge sweep angle of 0.4. The resiilts of the computations indicate that; 

1. For either the thickness or the volume condition; 

(a) The change in plan form resulted in only a sli^t change 
in the optimum section for either subsonic or supersonic leading 
edges . 

(b) The optimum airfoil section for supersonic leading edges 
is very nearly the two-dimensional optimum the drag reduction 
from the two-dimensional optimum is less than 10 percent. 

2. For the thickness condition: 

(a) Drag reductions of about 16 percent from double-wedge 
sections resulted for two examples having subsonic leading edges 
and supersonic maximum-thickness lines located at 50 percent of 
the chord. 
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(b) The optimum airfoil sections obtained for subsonic leading 
edges with subsonic maximum-thickness lines violate the assumptions 
of linear theory. Flow separation over these airfoils would be a 
governing factor in the drag determination. 

5. For the volume condition; 

(a) Drag reductions of 25 percent and 33 percent from biconvex 
sections resulted for a delta wing and an arrow wing, respectively, 
for a 60° swept leading edge at a Mach number of !.!<•. When com- 
pared with NACA 65A-series sections for the same conditions, the 
drag reduction is about 10 percent and 15 percent, respectively. 

(b) The optimum airfoil section for a very slender delta wing 
obtained by this analysis approaches the slender-body optimum 
section. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February I7, 195^* 
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TABLE I.- DELTA-WING AIRFOIL ORDINATES z/t AND DRAG PARAMETERS 

[» - I - 0-3 



2.357 


2.0677 


1.8412 


1.6745 

























CT\ 


XABIE rci— nBUTA-WIHG AJRKOIL (BDUMJEES z/t AHD IHAG PARAMEHEERS 

■= 0 } H “ 2 ^ 


(a) ThiolmesG condition (t) Aroa condition 



Ordinate^ z/t 



OirUnatflj z/t 

X 

= 

0.3 

”1 = 

1.0 

- 1.3 

= 

1.7 


X 

nj_ = 0.5 

nj_ ” 1.0 

Hi - 1.3 

ni “ 1.7 

n-i - 2,5 


5 - 0.2 

5-0.5 

1 - 0.2 

1-0,3 

0.2 

i - 0-2 

1-0.3 



K - 0.66312 

K - 0.64093 

K - 0.72793 

K - 0.72200 

E= 0.65085 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

.05 

.12232 

.04730 

.09318 

.03178 

.26726 

.17084 

.i&rr 5 


.05 

.08989 

.05477 

.26700 

.53747 

.35186 

.10 

. 2 U£ 1»-7 

.09527 

.21704 

.07379 

,56301 

.24078 

.25347 


.10 

.17165 

.12292 

.33477 

.41427 

.42793 

.15 

.57238 

.;iA386 

•35377 

.11981 

.43763 

.29181 

.31108 


.15 

.24496 

.19158 

.38351 

.46648 

.47^ 

.so 

.50000 

.19305 

.50000 

.16856 

.50000 

.50000 

.35685 


.20 

.30932 

.25707 

.41485 

.49090 

.49771 

.25 

.1)-6986 

. 843 ^ 

.40909 

.21948 

.46787 

.29353 

.39527 


.25 

.56507 

.51722 

.44063 

.50000 

.50000 

.30 

■J )-3951 

.29317 

.47172 

.27227 

.44557 

.26651 

.42737 


.50 

.41137 

.37051 

.46392 

.49903 

.49020 

.55 

. 1 ) 069 T 

.54407 

.44971 

.3^ 

,42145 

.^900 

.45353 


.35 

.44820 

.41577 

.48221 

.49191 

.47213 

. 1*0 

.37826 

.39551 

- 424 ® 


.39516 

.21935 

.47170 


.40 

. 4-^36 

.45^ 

.49448 

.48194 

. 44^3 

.Jv 5 

. 3‘^759 

.44749 

.39603 

.44063 

.36709 

.20106 

.48374 


■45 

.492^ 

.478a 

.50000 

,47126 

.42154 

.50 

.31637 

.50QOO 

■36571 

. 333 & 

.50007 

.50000 

.35756 

.18297 

.50000 


.50 

.50000 

. 4^77 

.49818 

.45789 

.44018 

.41732 

. 3 Sl^ 

.5921s 

■55 

.60 

.28^ 

.25393 

.43055 

.40091 

,45465 

.40778 

.30684 

.27512 

•i 6 te 7 

. 146 ^ 

.46322 

.42034 


.55 

.60 

.49718 

.48410 

.50000 

.49383 

II 

.36229 

.35311 

.65 

.22235 

.53117 

.26331 

.35^ 

.^259 

.12887 

.37368 


.65 

.46065 

,47587 

.44412 

.30380 

.70 

•15 

.19102 

.30151 

.22953 

.31047 

.20935 

.11071 

.32440 


.70 

.42675 

.44580 

.40872 

■35584 

.27^ 

. 1^0 

.25133 

.19266 

.26057 

.17552 

.09247 

.27318 


■P 

.58227 

.40535 

.56419 

• 31^5 

.23739 

.80 

.12769 

.20125 

.15544 

.20949 

.14119 

.07413 

.22047 


.80 

.32720 

.34852 

. 3 U «5 

.26419 

. 19 fe 

.85 


.15107 

.11736 

.15795 

.10641 

.05571 

.16660 


.85 

.26147 

.28053 

.24707 

.20691 

.153* 

.90 

, 0 ^ 1(00 

.10080 

.07872 

.10578 

.07126 

.C ^721 

.11179 


.90 

.18504 

.19994 

.17426 

. 14 ^ 

.10055 

.95 

.032014^ 

.05044 

.03958 

.05311 

.03577 

.01864 

.05621 


■95 

.09789 

.10642 

.09190 

.07685 

.05659 

1.00 

0 

0 

0 

0 

0 

0 

0 


1.00 

0 

0 

0 

0 

0 

PCd 

6.9175 

4.2211 

■Hi 

4.8363 

4.7317 

1.6412 

5.5906 


PCu 

t 2 

3.6416 

6.3176 

5-1339 

3.796B 

2.1993 
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TABIB 17 .- ARRCW-WIMG AIRFOIL OiEffilSAIES z/t AHD BRAG PARAMBimS 





0 ,l^} IT 



(a) Thdckneea condition (t) Area condition 




Ordinate j 

z/t 




Ordinate, z/t 

X 

“1 = 

0.5 

ni - 1.5 

ni « 1.7 

ni - 2.5 


X 

njL = 0.3 

“1 “ 1-3 

ni - 1.7 

Hi - 2.5 


1 = 0.2 

^=■0.5 

S = 0.5 

5 -O.S 

|o0.2 



K - 0.66423 

K - 0.73780 

K » 0.63403 

K = 0.60363 

0 

.05 

.10 
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TABLE V.- DELTA-WIKG AIREOIL ORDINATES BASED ON 
SLENDER-BODY THEORY 


X 

z/t~ 

0 

0 

.05 

.36060 

.10 

.45403 
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.25 

.48760 

.30 
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•95 
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1.00 
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Figure 2.- Exploded view of coefficient matrix of equations ( 9 ) to (l2) 


P\-g 







Figure 3.- Effect of number of chord, divtsions N on the drag and airfoil 
shapes of optimum delta wings. nSiickness condition; n]_ = 1 . 7 } ni/N = 0.20. 
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